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We study the Hubbard model with parameters relevant to the cuprates, using variational Monte Carlo with 
projected d-wave states. For doping < x < 0.35 we obtain a superconductor whose order parameter tracks 
the observed nonmonotonic T c (x). The variational parameter A va r(aO scales with the (tt, 0) "hump" and 
T* seen in photoemission. Projection leads to incoherence in the spectral function and from the singular 
behavior of its moments we obtain the nodal quasiparticle weight Z ~ x though the Fermi velocity remains 
finite as x — ► 0. The Drude weight D\ ovl and superfluid density are consistent with experiment and Di ow ~ Z. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.25.-q, 74.72.-h 
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Strong correlations are essential to understand <i-wave 
high temperature superconductivity in doped Mott in- 
sulators | p. The no-double occupancy constraint aris- 
ing from strong correlations has been treated within two 
complementary approaches. Within the gauge theory ap- 
proach P), which is valid at all temperatures, the con- 
straint necessitates the inclusion of strong gauge fluctua- 
tions. Alternatively, the constraint can be implemented 
exactly at T = using the variational Monte Carlo 
(VMC) method. Previous variational studies J3|-|| have 
focussed primarily on the energetics of competing states. 

In this letter, we revisit projected wavefunctions of the 
form proposed by Anderson in 1987 JjJ]. We compute 
physically interesting correlations using VMC and show 
that projection leads to loss of coherence. We obtain 
information about low energy excitations from the sin- 
gular behavior of moments of the occupied spectral func- 
tion. Remarkably, our results for various observables are 
in semi-quantitative agreement with experiments on the 
cuprates. We also make qualitative predictions for the 
doping (x) dependence of correlation functions in pro- 
jected states (for x <C 1) from general arguments which 
are largely independent of the detailed form of the wave- 
function and the Hamiltonian. 

We use the Hubbard Hamiltonian Ti = K, + Hint- 
The kinetic energy K. = J^k.a e ( k ) c Lr c ko- witn e ( k ) = 
— 2t (cos k x + cos k y ) + At' cos k x cos k y the dispersion on 
a square lattice with nearest (t) and next-near (t') hop- 
ping. The on-site repulsion is Hint = U Y^ r n-f(r)nj,(r) 
with n CT (r) = c£.(r)cv(r). We work in the strongly cor- 
related regime where J — 4t 2 /U, t' < t <C U near half 
filling: n = 1 — x with the hole doping x <C 1. We choose 
t' = t/4, t = 300meV and U = 12t, so that J = lOOmeV 
consistent with neutron data on cuprates. 

We describe the ground state by the wavefunction 



|*o> = exp(i5)P|*BCs). 
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Here |* BC s) = (E k ^( k ) c k T c -kiJ ' 1°) is the N ~ 
electron d-wave BCS wave function with y(k) = v^/uk 
A k /Kk 



/i V ar and A var determine <p(k) through £k = e(k) — /i V ar 
and Ak = A var (cosfca; — cosk y ) /2. The projector V = 
Y[ r (1 — n|(r)n^(r)) in Eq. ([!]) eliminates all configura- 
tions with double occupancy, as appropriate for U — > oo. 

The unitary transformation exp(iS) includes double 
occupancy perturbatively in t/U H|. The transfor- 
mation exp(— iS)HeyLj>(iS) is well known to lead to 
the tJ Hamiltonian M. We emphasize that our ap- 
proach effectively transforms all operators, and not only 
U. Using (* |Q|*o> = (*bcs|^Q^|*bcs) with Q = 
exp(—iS)Q exp(iS), it follows that incorporating exp(iS) 
in the wavefunction (1) is equivalent to transforming 
Q—> Q and working with fully projected states. 

Using standard VMC techniques || we compute 
equal-time correlators in the state |^o)- The two 
variational parameters are determined by minimizing 
(*o|W|*o)/(*o|*o) at each x. The doping dependence of 
the resulting A var (:r) is shown in Fig. 1(a). Varying input 
parameters in TL we find that the scale for A var is mainly 
determined by J = 4i 2 /U. We show below that A var is 
not proportional to the SC order parameter, in contrast 
to BCS theory, and also argue that it is not equal to the 
spectral gap. On the other hand, we find that the optimal 
Mvar(a^) ~ Mo (a;) j the chemical potential of noninteracting 
electrons described by /C. However, ^ var is quite differ- 
ent jj| from the chemical potential ji = d(H)/dN, where 
(• • •) denotes the expectation value in the optimal, nor- 
malized ground state. 

Phase Diagram: We now show that the wavefunc- 
tion (Q) is able to describe three phases: a resonating 
valence bond (RVB) insulator, a d-wave SC and a Fermi 
liquid metal. To establish the T = phase diagram we 
first study off-diagonal long range order (ODLRO) us- 



A£]. The two variational parameters 



ing F a< p(r — r') = (cl(r)cj(r + d)cx(r')c T (r' + $)), where 
a,j3 = x, y. We find that F a ,p — » ±$ 2 for large |r — r'|, 
with + (— ) sign obtained for a || (_L) to [3, indicating 
d-wave SC. As seen from Fig. 1(b) the order parameter 
<$>{x) is not proportional to A var (a;), and is nonmonotonic. 
<I> vanishes linearly in x as x — ► as first noted in ref . [EJ , 
even though A var ^ 0. We argue |l(| that $ ~ x is a 
general property of projected SC wavefunctions. The lo- 
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cal fixed number constraint imposed by V at x — leads 
to large quantum phase fluctuations that destroy SC or- 
der. The non-SC state at x — is an insulator with a 
vanishing Drudc weight, as shown below. The system is 
a SC in the doping range < x < x c ~ 0.35 with <E> ^ 0. 
For x > x c , $ = and the wavefunction ^>q for A var = 
describes a normal Fermi liquid. In the remainder of this 
paper we will study < x < x c . 
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FIG. 1. (a): The variational parameter A var (filled squares) 
and the (n, 0) hump scale (open triangles) in ARPES Jill] 
versus doping, (b): Doping dependence of the d- wave SC 
order parameter $. Solid lines in (a) and (b) are guides to 
the eye. (c): The coherence length £ sc > max(£ pa i r , 1/ y/x). 



Coherence Lengths: We must carefully distinguish 
between various 'coherence lengths', which are the same 
in BCS theory, but are very different in strongly corre- 
lated SC's. The internal pair wavefunction ip(k) defines 



cal correlations can be expressed as equal-time correla- 
tors, calculable in our formalism. Specifically, we look 
at M^(k) = du)(j e f(u))A(k, uj), where A(k,u>) is the 
one-particle spectral function and, at T = 0, f(u>) — 
@(—u>). We calculate Mn(k) = n(k) and the first mo- 
ment jy| Mi(k) = (cj^JW, Cko-]) to obtain important in- 
formation about the spectral function. 

Nodal Quasiparticles: Fig. 2(c) shows that along 
(0, 0) to (-7T, 7r) n(k) has a jump discontinuity. This 
implies the existence of gapless quasiparticles (QPs) 
observed by angle resolved photoemission spectroscopy 
(ARPES) |Q. The spectral function along the diagonal 
thus has the low energy form: A(k,ui) = Z8(u — £&) + 
Ai nc , where = vp(k — hp) is the QP dispersion and 
Ai nc the smooth incoherent part. We estimate the nodal 
kp(x) from the location of the discontinuity and the QP 
weight Z from its magnitude. While kp(x) has weak dop- 
ing dependence, Z{x) is shown in Fig. 2(d), with Z ~ x 
as the insulator is approached [jl5| . 

Projection leads to a suppression of Z from unity with 
the incoherent weight (1 — Z) spread out to high energies. 
We infer large incoherent linewidths as follows, (a) At the 
"band bottom" n(k = (0,0)) ~ 0.85 (for x = 0.18) im- 
plying that 15% of the spectral weight must have spilled 
over to u> > 0. (b) Even at kp, the first moment Mi 
lies significantly below u) = (defined by the chemical 
potential fj, = d(H)/dN); see Fig. 3(a). 

The moments are dominated by the high energy inco- 
herent part of .A(k, ui), but their singular behavior is de- 
termined by the gapless coherent QPs. Specifically, along 
the zone diagonal Mi (k) = Z£k®(—£,k) + smooth part. 
Thus its slope dMi(k) / dk has a discontinuity of Zvp at 
fcp, as seen in Fig. 3(a), and may be used to estimate 
]l6[ the nodal Fermi velocity vp. As seen from Fig. 3(b), 



_n / » i -n • ii i vp(x) is reduced from its bare value d p and is weakly 

Vp A var , where is the bare average , . , , . ., , .^.-r-.™ . 



a pair-size £ pair 
Fermi velocity. Projection does not affect the pair-size 
much, and £ pa i r remains finite at x — 0, where it defines 
the range of singlet bonds in the RVB insulator [D. 

A second important length scale is the inter-hole spac- 
ing 1 / y/x. At shorter distances there are no holes and the 
system effectively looks like the x — insulator. Thus 
the superconducting coherence length £ sc must necessar- 
ily satisfy £ sc > max (£ pa i r , 1 / y/x) . This bound implies 
that £ sc must diverge both in the insulating limit x — > 
(see also Refs. 0) and the metallic limit x — * x~ , but is 
small at optimality; see Fig. 1(c). This non-monotonic 
behavior of £, sc (x) should be checked experimentally. 

Momentum Distribution: From gray-scale plots 
like Fig. 2 we see considerable structure in the momentum 
distribution n(k) over the entire range < x < x c . One 
cannot define a Fermi surface (FS) since the system is a 
SC (or an insulator). Nevertheless, for all x, the k-space 
loci (i) on which n(k) = 1/2 and (ii) on which |VkH(k)| 
is maximum, are both quite similar to the corresponding 
noninteracting FS's. For t' = t/4, the FS is hole-like for 
x < 0.22, and electron-like for overdoping Jl2[ . 

We next exploit the fact that moments of dynami- 



doping dependent, consistent with the ARPES estimate 
of vp w 1.5eV-A in Bi 2 Sr 2 CaCu20 8+5 (BSCCO). 

As x -> 0, Z = {l-dZ'/duj}- 1 ~ x, while v F (x)/v F = 
Z\l + (v F )~ l dYi 1 / dk] is weakly x-dependent. This implies 
a 1/x divergence in the k-dependence of the self energy S 
on the zone diagonal, which could be tested by ARPES. 

Within slave boson mean field theory (MFT) §, we 
find H Z sh =x in Fig.2(d) and vf(x) shown in Fig.3(b), 
both systematically smaller than the corresponding VMC 
results. Thus the holons and spinons must be partially 
bound by gauge fluctuations beyond MFT. 

Moments Along (7r,0) — > (tt,tt): The moments n(k) 
and Mi(k) near k = (it, 0) are not sufficient to estimate 
the SC gap, however they give insight into the nature 
of the spectral function. As seen from Fig. 4(a) n(k) is 
much broader than that for the unprojected I^bcs)- For 
k's near (n, 7r), which correspond to high energy, unoc- 
cupied (u> > 0) states in I^bcs) we see a significant build 
up of spectral weight transferred from uj < 0. Corre- 
spondingly, we see loss of spectral weight near (ir, 0). We 
thus infer large linewidths at all k's, also seen from the 
large |Mi(k)| (in Fig. 4(b)) relative to the BCS result. 
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The transfer of weight over large energies inferred 
above suggests that projection pushes spectral weight to 
\uj\ < A var , the gap in unprojected BCS theory. We thus 
expect that the true spectral gap -E gap < A var and that 
the coherent QP near (tt, 0) must then reside at E gap in 
order to be stable against decay into incoherent excita- 
tions. It is then plausible that A var , the large gap before 
projection, is related to an incoherent feature in A(k, uj) 
near (tt, 0) after projection p7| . Indeed, comparing A var 
with the (tt, 0) "hump" feature seen in ARPES 11 , we 
find good agreement in the magnitude as well as doping 
dependence; see Fig. 1(a). 

Optical spectral weight: The optical conductivity 



sum rule states that J °° dujRea(uj) = nJ2i 



L (k)n(k) 



= nD tot /2 where m _1 (k) = (d 2 e(k)/d\i x d\i x ) is the 
noninteracting mass tensor (we set fi = c = e = 1). The 
total optical spectral weight D to t(x) is plotted in Fig. 5(a) 
and seen to be non-zero even at x = 0, since the infinite 
cutoff in the integral above includes contributions from 
the "upper Hubbard band" . 

A physically more interesting quantity is the low fre- 
quency optical weight, or Drude weight, D\ ow where 
the upper cutoff extends above the scale of t and J, 
but is much smaller than U. This is conveniently de- 
fined Q by the response to an external vector potential: 
Aow = d 2 (TiA)A/dA 2 . Here the subscript on TL denotes 
that A enters the kinetic energy via a Peierls minimal 
coupling, and that on the expectation value denotes the 
corresponding modification of exp(iS) in Eq. (^). The 
Drude weight D\ ow (x) 1 plotted in Fig. 5(a), vanishes lin- 
early as a; — > 0, which can be argued to be a general 
property of projected states ||. This also proves, follow- 
ing Ref. Jig] , that \^o) describes an insulator at x = 0. 
Both the magnitude and doping dependence of D\ ow (x) 
are consistent with optical data on the cuprates jl9) . 

Motivated by our results on the nodal Z(x) and 
D\ ow (x), we make a parametric plot of these two quanti- 
ties in Fig. 5(b) and find that D\ ow ~ Z over the entire 
doping range, a prediction which can be checked by com- 
paring optics and ARPES ||(J on the cuprates. 

Super-fluid Density: The Kubo formula for the su- 
perfluid stiffness D s can be written as D s = D\ ow — A T , 
where D\ ow is the diamagnetic response, and A T is the 
transverse current-current correlator. Using the spectral 

> 
First, 
22). Sec- 



representation for A T it is easy to see that [^lj A 1 



< Di ow . Two conclusions follow 
0, consistent with experiments 



implying D s 
D s — > as x 

ond, we obtain a lower bound on the penetration depth 
A L defined by A^ 2 = 4:ire 2 D s /h 2 c 2 d c . Using the calcu- 
lated -Di ow ~ 90meV at optimality and a mean inter- 
layer spacing d c = 7. 5 A (appropriate to BSCCO), we 
find A L > 1350 A, consistent with experiment p3[ . 

Conclusions: We have shown, within our variational 
scheme, the T — evolution of the system from an un- 
doped insulator to a d-wave SC to a Fermi liquid as a 
function of x. The SC order parameter <S>(x) is non- 
monotonic with a maximum at optimal doping x ~ 0.2, 
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FIG. 2. (a) and (b): Gray-scale plots of n(k) (black = 1, 
white = 0) centered at k = (0, 0) for x = 0.05 and x — 0.18 re- 
spectively on a 19 x 19 + 1 lattice, showing very little doping 
dependence of the large "Fermi surface", (c): n(k) plotted 
along the diagonal direction (indicated as a in Panel (b)), 
showing the jump at Icf which implies a gapless nodal quasi- 
particle of weight Z. (d): Nodal quasiparticle weight Z(x), 
compared with the slave boson mean field Z Bb (x) (dashed 



line) . 
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FIG. 3. (a): The moment Mi(k) along the zone diagonal, 
with smooth fits for k < Iqf and k > I(f, showing a disconti- 
nuity of Zvf in its slope at kF- (b): Doping dependence of 
the nodal quasiparticle velocity obtained from Mi(k). Error 
bars come from fits to Mi(k) and errors in Z. Also shown are 
the bare nodal velocity v%, the slave boson mean field fj^a:) 
(dashed line), and the ARPES estimate ?4° xpt) @. 



suggestive of the experimental trend in T c (x) from under- 
to over-doping. As x — > D s vanishes, while the spectral 
gap, expected to scale with A var , remains finite. Thus the 
underdoped state, with strong pairing and weak phase co- 
herence, should lead to pseudogap behavior in the tem- 
perature range between T c (which scales like $) and T* 
(which scales like A var ). 
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FIG. 4. (a): n(k) plotted along the (n, 0)-(7T, 7r) direction 
(indicated as /3 in Fig. 2(b)) and compared with the BCS 
result, (b): The moment Mi(k) plotted along the (tt, 0)-(tt, 7r) 
direction compared with the BCS values. 
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FIG. 5. (a): Doping dependence of the total (Dtot) and 
low energy (Di ow ) optical spectral weights (b): The optical 
spectral weight Di ow versus the nodal quasiparticle weight Z. 
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